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Abstract 

Oh! 

(— I ' We compute by Bethe Ansatz both bulk and boundary hole scattering matrices for 

'"^ ' (1) 

the critical A\/_i quantum spin chain. The bulk S matrix coincides with the soliton 

K^ \ S matrix for the Aj^_-^ Toda field theory with imaginary coupling. We verify our 

^ ■ result for the boundary 5 matrix using a generalization of the Ghoshal-Zamolodchikov 



boundary crossing relation. 



1 Introduction 

Magnetic chains associated with affine Lie algebras [|l|, constitute large classes of inte- 
grable models. These models hold both mathematical and physical interest. Indeed, the 
models can be shown to have certain quantum-algebra symmetries which are characterized 
by a parameter q. A variety of powerful analytical techniques can be brought to bear on 
these models, yielding in the thermodynamic limit exact properties such as excitation spec- 
tra, scattering matrices, and correlation functions. These models, which can be regarded 
as discretized integrable quantum field theories, have numerous applications to statistical 
mechanics and condensed matter physics. 

In a previous paper |^ we began investigating the simplest such class of models, namely, 
the ^V-i ^P^^ chain. This class includes the well-known XXZ chain as the special case 
A/" = 2. There we considered the noncritical regime (real g), for which there is a mass gap, 
and for which the S matrices are characteristic of g-deformed quantum field theory. 

In the present work we continue our investigation of the ^V-i ^P^^ chain, but consider 
instead the critical regime (|g| = 1). We compute by Bethe Ansatz both bulk and boundary 
hole scattering matrices. Our result for the bulk S matrix coincides with the soliton S 
matrix for the ^V-i Toda field theory with imaginary coupling [§]-0. We also generalize 
the Ghoshal-Zamolodchikov boundary crossing relation to the ^V-i ^^'^^ (^°^ which the 
bulk S matrix does not have crossing symmetry), and use it to help verify our result for the 
boundary S matrix. For A/" = 2, the boundary S matrix reduces to the well-known result 
for the critical XXZ open spin chain [0-|0; while in the isotropic limit, it coincides with the 



result found in 10 



2 The model and its symmetries 

We begin by briefly reviewing the construction of the model, both closed (periodic boundary 
conditions) and open with diagonal boundary fields. Q For an introduction to the Yang- 



Baxter equation, integrable spin chains, and the Bethe Ansatz, see e.g. [O], ||T2 



2.1 Closed chain 

The main ingredient for constructing an integrable closed spin chain is a solution R of the 
Yang-Baxter equation 

i?12(A) i?l3(A + A') i?23(A') = i?23(A') i?13(A + A') i?i2(A) . (1) 

We consider here the ^V-i -^ niatrix [^ 



R{\) = aY, tjj ® e,-, + h J2 ejj ® Ckk + cJ2 ^jk ® e^, e"'^^ ^*3"(^-^) , (2) 

^Our notations closely follow those of Ref. Q, in which the noncritical regime is studied, with the 
anisotropy parameter rj replaced by i/i. 



where the indices j , k take values from 1 to N", and Cjk are elementary N" x Af matrices with 
matrix elements (cja:)^^ = SjaSkb- Moreover, 

a = sinh (/i(A + i)) , b = sinh(/iA) , c = sinh(i/i) , (3) 

where the anisotropy parameter /i is real. The R matrix becomes SU{A/) invariant for /i ^ 0, 
and is proportional to the permutation matrix for A = 0. 

The corresponding closed-chain transfer matrix t(A) is given by 

t(A)=troTo(A), (4) 

where trg denotes trace over the "auxiliary space" 0, and Tq{X) is the monodromy matrix 

To(A)=i?oiv(A)---/?oi(A). (5) 

We restrict the number of spins, A^, to be an even number. 

The closed-chain Hamiltonian 7i is given by the logarithmic derivative of the transfer 
matrix at A = 0, 



Hr^— logt(A) ^^^ = -(5: /i„,„+i + /i^,i) , (6) 

n=l 



d\ 

where the two-site Hamiltonian h is given by 



Here we restrict yU to lie in the range /x G (0 , tt). As for the ^V-i Toda field theory with 
imaginary coupling, the Hamiltonian for A/" > 2 is not Hermitian. Nevertheless, we find in 
our study of low-lying states (see Sec. 3) only real energy eigenvalues. 

The transfer matrix has an exact U{1) symmetry. Indeed, let 

s^^^ = ekk-ek+i,k+u A; = 1,... ,7V- 1, (8) 

be the Cartan generators in the defining representation of SU{J\f)] and let s^^^ be the corre- 
sponding generators at site n, 

i 

s^f^ = 1®...®!® sW ®l®...(g)l, n = l,...,N. (9) 

We denote by S^''^ the corresponding "total" generators acting on the full space of states 

N 

s^'^ = J2sl,'K k = i,...,M-i. (10) 

n=l 



^The Hamiltonian with /i e (tt , 27r) is related to the parity-transformed Hamiltonian with ^ G (0,7r). 
Explicitly, 7i(/x) = H7i(/i')H, where ii' — 2tt — /i, and H is the parity [^ operator. An alternative approach 
(followed, e.g., in |8),[g[) is to restrict /i even further to the range /^ e (0, f) and to consider both signs of 
the Hamiltonian. 



The transfer matrix has the symmetry 

't{X),S^^A =0, 



k= 1,... ,A/'- 1. 



;ii) 



The eigenstates and eigenvalues of the transfer matrix have been determined |114| , [p!5[l , |12 



by the nested algebraic Bethe Ansatz method. Indeed, the states are constructed using 
certain creation operators depending on the solutions {A^^} of the Bethe Ansatz equations 



M(J-i) 



M(J) 



1 = n e.^{X^£-Xt''■,^.)l[e2{X^<^^-V■^^^) IT e_i(A(/-) - A^^+^^ /i) 



(i). 



M(J+i) 



.(i+1). 



(3=1 



0=1 



/3=1 



a 



i,...,M(^), J = l,...,^f-l, 



where 



e„(A;yu) 



sinh/i(A + f) 
sinh/x(A-f) ' 



(12) 



(13) 



and M^^-* = N , M^-^'^^'^ = , A^''-' = A^-* = . The corresponding energy, momentum, and 
gik) eigenvalues are given by 



mW 



E 



— sm' 



>E 



(l)^ 



1 cosh(2/iAa ) — COS/i 



mw 



P = -^logei(A«;/i) (mod27r) 



«=i 



g{k) ^ Af (^"^) + M(^+^) - 2M^'''> 



(14) 

(15) 
(16) 



2.2 Open chain 

In addition to an R matrix, the construction of an integrable open spin chain requires also 



a solution K of the boundary Yang-Baxter equation 16 



RuiXi - X2)K,{X^)R2l{X^ + X2)K2{X2) = K,{X2)Ru{X, + X2)K^{X,)R2l{X^ - X^) . (17) 
We consider here the diagonal ^V-i ^ niatrix |[TP| 



3=1 j=l+l 



where 



a = sinh {^{iC, — A)) , /5 = sinh {^{iC, + A)) , 



(19) 



where ^ is a parameter and I E {1 , . . . ,Af — 1}. 

The corresponding open-chain transfer matrix t(/)(A ,,^_ ,(,+) is given by p!7|-||T9|,HI 

t(0(A , e- , e+) = tro Mo K^i) o(-A - P , e+ - y ) To(A) ir(,) o(A , ^.) fo(A) , (20) 

where the monodromy matrix To(A) is given by (j^), and Tq{X) is given by 

fo{X)=Rio{X)---Rm{X). (21) 

Moreover, M is the matrix in the crossing-unitarity relation 

Ri2{-P - X)''MiRu{-p + XY'M{^ oc 1 , (22) 

and is given by the J\f x J\f matrix Mj^ = Sj^e^^'^ '^^^^\ and p = iJ\f/2. Indeed, it can be 
shown that this transfer matrix has the commutativity property 



_t(o(A,e-,e+),t(/)(A',e-,e+) 

and its derivative at A = gives the Hamiltonian, 



0, 



'H ~ Tx'^^')'-'^'^"'^"^ 



N-l 



A=0 2 



J2 hn,n+l + C(e-)P(/) 1 - C(e+ -^f+l)Pi 



(/) TV 



(23) 



(24) 



n=l 



where the two-site Hamiltonian h is given by (^, and 



c(0 = -sin/i(i-cot(/iO) 






^(0 = H ^31 



Z^ ^ij 



(25) 



In particular, for M = 2 (and therefore / = 1), the Hamiltonian is given by [^ ,117 



-^ N-l 
^ = 7(11 («+i + <^n+i + cos/i«+i) - sin/icot(/i^_)(T^ + sin /i cot (/i(^+ - l))a^} .(26) 

^ n=l 

Evidently, the parameters ^-t correspond to boundary fields. 

An important feature [Q of this model is its exact quantum algebra symmetry Uq {SU{1)) 
X Ug {SU{Af — I)) X f/(l). Indeed, let J^'^^^ {k = 1 , . . . , A/"— 1) be raising/ lowering operators 
of the quantum algebra Ug {SU{J\f)) which act on the full space of states, and which obey 
the commutation relations 



j+(fc) , j~U) 



Sk,j 



■^(fc) 



Sik) ^ j±0-)J = ± ^26kj - 6k-i, - 4+1,,) J^^'^ , (27) 

where [x]^ = {q^ — q~^)/{q — q^^), and S^'^^ is given by Eq. (p^O]). For a given value of /, 

t(0(A,e-,e+),5('^] = 0, k = l,...,Af-l, 

= 0, A;^/, (28) 

4 



t(/)(A,^,e+),J^(') 



where q = e *^. Had we taken the K matrix to be the identity matrix K{X) = 1 instead of 



[^), then the corresponding transfer matrix would have the full Uq {SU{N')) symmetry [|T8 
The model also has P] a "duality" symmetry which relates / ^^ JV — I, 



W t^i) (A , e- , i+) U-' ex ^M-i) (A , -i- , -i+ + M) , 



(29) 



i(i) 



where lA = 1. This is a remnant of the cyclic (Zj^) symmetry |T^, [^ of the Aj^'_i R 
matrix. 



The Bethe Ansatz equations are [|19 

X n e_i(A(^-)-Aj/--^);^)e_i(A(^VA(^--^);/i) ne2(AL^^)-A(^-);/i)e2(A(^VA(^^^ 



mO+1) 

X n e_i(A(f)-A(^--^^);/.)e_,(A(^-) + A(^--^^);/z) 

/3=1 



/3=1 



a 



1,...,M(^-), j = l,...,Ar-l. 



(30) 



The energy is given by ([T^) (plus terms that are independent of {A^^^}), and the S^''^ eigen- 
values are again given by Eq. (|1^). 



3 Bulk S matrix 



In order to investigate bulk properties, we consider the periodic chain with Bethe Ansatz 
Eqs. (|T^. We assume that the ground state is the Bethe Ansatz state with no holes, i.e., 
with Af — 1 filled real Fermi seas. Holes in these seas correspond to "solitons" . Indeed, the 
Bethe Ansatz state with one hole in the j*'' sea is a particle-like excited state which belongs 
to the fundamental representation [j] of Uq {SU{Af)), corresponding to a Young tableau with 
a single column of j boxes. Such an excitation with rapidity A^-'-' has energy ^^™^ g(j)(A(-^)) 

and momentum ^'^•'^(A'--''*), where 



s^'^iX) 



and the momentum satisfies 



sin(^(Ar-j)) 



Arcos(^(AA-j)) + cosh(^)' 



^A,0)(,).,0)(,) 



(31) 



(32) 



Note that s^-'^(A) has the periodicity A — > A + ij\f, which suggests that the physical strip is 



< 3m A < 



Af 



Exact bulk scattering matrices can be computed using a generalization of the method of 
Korepin [^ and Andrei-Destri ||23|. We define the two-particle S matrix Rbl^t^l for particles 



of type [j] and [k] by the momentum quantization condition 

^^^p(^Hm)N^mm _ ij \~x^j) ^ 3,W) = o . (33) 

For the scalar factor, we obtain (cf. [|10]) 



exp 2 r — sinh (z^(A(^) - A^^'))) {5,k - RA^)) , (34) 



00 



where a^^\\) is the density of Bethe Ansatz roots {A^^} for the state with holes of rapidities 
A^-'^ and A*^'^-' in the j*'' and /c*^ seas, respectively; and (see, e.g., |jl2[) 

sinh(f^) sinh(^) sinh ((TV- j>)f ) 
^^ sinh ((z/- l)f) sinh (^) sinh (f) ' 

where u = -, j> = max(j ,j') and j< = min(j , j'). Note that yU G (0 , tt) implies z/ > 1. 

Let us focus on the particular case of two holes of type [1]. The scalar factor is given by 

rW^W -exvl2r- sinh Uu;~X] ^^^^ (f ) sinh ((z. - AT)^^^ 
' - '"P f io .; '"^ i^^^^ sinh ((. - 1)1) sinh (^'^ ^ ' ^''^ 

where A = A^ — A2 • Moreover, let us consider the following Bethe Ansatz states: 

(a) two holes and one 2-string (4-string) of positive parity in the first sea if z/ > | ( z/ < |, 
respectively) 

(b) two holes and one 1-string (3-string) of negative parity in the first sea if z/ > 2 ( z/ < 2, 
respectively) 

The corresponding S matrix elements are given by 

It still remains an open problem to determine all N"^ two-particle states. However, let us 
assume that the scattering matrix R[^]®[^1 has the same structure as the R matrix (0); i.e., 

R(A) W« W = a Y. (^n ® e,, + b J2 ^n ® e^^ + c ^ e,k ® e^, e"^'^ ^^^'^^^-'^^ , (38) 



■^Up to signs, which are not easy to obtain from the Bethe Ansatz, but which can be fixed by requiring 
that for A = the amphtudes be equal to +1 and —1 respectively; i.e., that the matrix R(0) (see below) be 
equal to the permutation matrix. 



where a , b , c , x are to be determined. This matrix has the three distinct eigenvalues a , b + 
c , b — c, which we identify with the three amphtudes Rg , R(a) , R(b) , respectively. Note 
that the eigenvalues of R[^]®[^] are independent of x. That is, the value of x cannot be 
determined by considering Bethe Ansatz states with two holes. We fix the value of x (up to 
a sign) by requiring that R[^l®[^] obey the Yang-Baxter Eq. (|1]).0 In this way, we obtain 

a = sin ( -{1 + iX))r{X) , b = — sin ( -jr(A) , c = sin ( -jr(A) , x= r , (39) 



where 



r(A) = ^ ^. (40) 

sin(^(l+-^ 



For A/" = 2, this S matrix agrees with the known result for the XXZ chain [^, |]TB[. Moreover, 
the scalar factor (|3^) coincides with the expression given in Eq. (3.29) of Ref. [^ for the 
a^^ scalar factor in afiine Toda field theory with imaginary coupling, provided we make the 
following identifications: 

n^Af-1, \^ -, /x^^— . (41) 

z/ — 1 u — I 

It has been observed |^-|§] that (after a certain gauge transformation) this S matrix 
has a Uq{SU{N')) symmetry with q = e'-^^ . On the other hand, from (p8|) we see that 
the finite-size transfer matrix has an "approximate" Uq{SU{N')) symmetry with q = e~ . 
Evidently, as a result of filling the Fermi seas and taking the thermodynamic limit, the value 
of q becomes renormalized. (See also e.g. 



4 Boundary S matrix 

We define the boundary S matrices K|^-j r .i for a particle of type [j] by the quantization 
condition [^ 

(e"^"'^'"'^^Kj,,K,),,-l)|A(^-)) = 0. (42) 

For simplicity, we consider only the case of a hole of type [1]. The quantum-algebra symmetry 
(p8|) implies that the boundary S matrices K|^> m, are diagonal A/" x A/" matrices of the same 
form as the K matrix (0), 

K^) [1] = "^)e^^^"' E e,, + P^,f^y''''' E ^n , (43) 



^Presumably, one can alternatively consider Bethe Ansatz states with three holes, and require factorization 
of the three-particle S matrix into a product of two-particle 5' matrices, as was done for the XXX model in 



where aj^x, /3Jn, y are to be determined. We have the relation 

a+ a^) ~ exp |^27riV j^' ' (a«(A) - 2s^'\\)) d\\ , (44) 

where anHX) is the density of Bethe Ansatz roots {A^^^} for the state with one hole of 
rapidity A^^^ in the first sea. From the Bethe Ansatz Eqs. (|30D , we find (cf. ||3|, [p^] ) 

a^)=zcosh(-^(A«-^(z.-2^))) A:(A«,e-), (45) 

where 

k(\,i) = ko{X)HX,0, (46) 

with 

( roodiu / ^, sinh ({M + l)f ) sinh ({y - M)'^^ 

UX) = exp 2 / — sinh(2z^A) ^^ -^ ^ '^ ^ 

[ Jo uj V y sinh (ATtu) sinh ((z/-l)f 

hCx,0 = ^ 



;cosh(^(A-i(i/-2e 

f ^ f^ du . ^ / . ~N sinh((Ar-/)cu) sinh ((z/ - 2^ - /)a;) 1 ,^^, 

X exp I -2 — smh {2tujX] ^ , ,/ \ x , ..r s ^ • (47) 

^\ io w ^ ^ sinh((z/-l)cj) sinh(Arw) J ^ ' 

Moreover, with the help of the "duality" symmetry (^9]) , we obtain 

/3^) = zcosh(-^(A« + ^(z.-2U)) A:(A('^e-)• (48) 

Our approach for computing a7i^ and (3n^ involves performing certain Fourier transformations, 
which leads to the following restriction of parameters: J\f < 2^_ + / < 2j/. In particular, this 
requires u > y- 

The parameter y can be determined using the boundary Yang-Baxter Eq. ([17D (in parallel 
with our analysis of the parameter x in the bulk 5* matrix), and we obtain the value 

(49) 



u-1 



Our expression for the boundary S matrix agrees with [0-0 f\ for A/" = 2, and with IJID 
in the isotropic limit z/ — i> oo. A further check on our result is performed in the following 
section. Finally, we remark that for Kt-, m, there are similar expressions involving C,+- 

'''The factors e^"^ in the boundary S matrices were unfortunately omitted in H. 
^In H we use a slightly different definition of the boundary parameters ^ip. 



5 Functional equations for boundary S matrix 

Thus far, we have pursued the direct Bethe-Ansatz calculation of S matrices. Alternatively, 
S matrices can be found (up to CDD ambiguities) by the "bootstrap" approach |2^ . Indeed, 
the bulk S matrix scalar factor (^) was first obtained for affine Toda theory in this way 

i-i- 

A natural question is whether the boundary S matrix can also be obtained by a boot- 
strap approach. For the case that the bulk 5* matrix has crossing symmetry, Ghoshal and 
Zamolodchikov ^ have formulated a boundary crossing relation, which - together with the 
boundary unitarity relation - determines the scalar factor of the boundary S matrix, up to 
a boundary CDD ambiguity. Unfortunately, the ^V-i bulk S matrix does not have crossing 
symmetry for A/" > 2, and hence, the Ghoshal-Zamolodchikov relation cannot be directly 
applied. 

In this section, we generalize the Ghoshal-Zamolodchikov boundary crossing relation to 
the AjiJ_i case, and we find that the boundary S matrix scalar factor P6|),(^) can indeed 
be obtained by this bootstrap approach. The basic observation is that, although the bulk S 
matrix lacks crossing symmetry, it does have crossing-unitarity symmetry; and one expects 
that a corresponding property should hold for the boundary S matrix. [] 

We begin by briefly reviewing the case that the bulk S matrix does have crossing sym- 
metry, and rewriting the Ghoshal-Zamolodchikov boundary crossing relation in matrix form. 
We assume that the bulk S matrix R(A) has PT symmetry, 

Ri2(A)*^*^ = R2i(A), (50) 

where R21 = 'Pi2Ri2'Pi2) T^ is the permutation matrix, and tj denotes transposition in the i*'^ 
space. We write the bulk crossing relation as 

Ri2(A)*^ = riRi2(p-A)Vi, (51) 

where the crossing matrix V satisfies V"^ = 1, and p is the crossing parameter. The unitarity 
relation is 

Ri2(A)R2i(-A) = l. (52) 

Combining the three relations (^)-(|5^), one obtains the crossing-unitarity relation 

Ru{p-~XY'M,Ru{p+~Xy'M^' = l, (53) 

where M = V^V. We observe that the Ghoshal-Zamolodchikov boundary crossing relation 
R] for the boundary S matrix K(A) is consistent with 



tri 



R2i(-2A)MiKi(p + A)l = (v^K^i-mY" , (54) 



^A different generalization of the Ghoshal-Zamolodchikov boundary crossing relation has recently been 
discussed in Ref. [E9| . 



where R12 = ^12^12, and hence R21 = ^u'Pu- (For aesthetic reasons, we have made a shift 
in the rapidity variable by p/2, but this is not necessary.) 

With the help of the boundary unitarity relation 



K(A)K(-A) = 1 



one can now obtain the relation 



tr 



13 



R2i(-2A)MiKi(p + A)R23(2A)M3K3(p - A) 



1. 



(55) 



(56) 



We shall see below that this relation, which is a boundary generalization of the bulk crossing- 
unitarity relation (|53|), is the desired result. 

Let us now turn to the ^V-i bulk S matrix R[i]®[i]. Although it does not have crossing 
symmetry, it does obey the crossing-unitarity relation (|5^) with Mjk = Sjke~^^ ~ "'^ >'i^^^ 
and p = iAf/2 (cf. Eq. (|2^) ). We observe that this relation together with the unitarity 
relation ( |52D determine the scalar factor of the bulk S matrix. Indeed, substituting the form 
(^), (^) into these relations, we obtain the following functional equations for r(A): 

,iAf r, ,iAf r, 1 

^(^ + A)r(— -A) 



sm 



_n_(f£ 
1/-1V 2 



r(A) r(-A) 



+ a))sin(^( 
1 



2 



iX) 



sin -^(1 + ^A) sin^;^(l-^A) 



(57) 



Solving these equations yields the result (^OD, (|36D, up to a CDD factor. 

Finally, we turn to the AjJ_i boundary S matrix K/^^ r-^,. Substituting the form (|i5|), (^), 
(^) into the boundary crossing-unitarity relation (|56|), we obtain a functional equation for 
the scalar factor k{\ , ^): 

H^ + A nn^-A ^i^te(i + 2^^))^i^te(i-2^^)) 

2 '^^ ^ 2 '^' sin(-^(Ar+2zA))sin(^(Ar-2a) 



-11 



X 



cosh (^(A - |(z/ - 2^ + TVT - 2/))) cosh (^(A + ^{u-2^ + Af- 21)) 



{51 



(up to a rapidity-independent phase, which can be absorbed by a redefinition of M). Similarly, 
the boundary unitarity relation (EH) implies 



1 



cosh (^(A - f (z/ - 20)) cosh (^(A + f (z/ - 2^ 
Solving these equations yields the result (ffSp, (^Tj) up to a CDD factor. 



(59) 



10 



6 Discussion 

The critical ^V-i quantum spin chain shares a number of features with ^V-i Toda quantum 
field theory with imaginary coupling. In particular, the two models have solitonic excita- 
tions with the same bulk 5* matrix. This suggests that the spin chain may be a type of 
discretization (see, e.g., [^,|^) of the affine Toda theory. This further suggests that it 



may be possible to find "soliton preserving" integrable boundary conditions for the affine 
Toda theory which lead to the boundary S matrix which we have given here. Searches for 
such integrable boundary conditions have not yet been successful (see, e.g., |^,|^), but 



progress in this direction has recently been reported |^. (Boundary S matrices with a 
structure similar to ( ^31) have recently been proposed for the 0{N) Gross- Neveu and nonlin- 



ear sigma models [^ and for the SO{N) principal chiral model [^.) As in the affine Toda 



theory, the spin chain has various bound states, on which we hope to report soon. 
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Note Added: After completion of this work, we became aware of Ref . |^ , which puts 
forward the more general conjecture that the light-cone continuum limit of the critical g spin 
chain is the g Toda theory with imaginary coupling, for any untwisted affine simply-laced 
Lie algebra g. The result (^) for the scalar factor of the bulk 5* matrix for the ^V-i chain is 
also obtained there. We are grateful to G. Takacs for bringing this reference to our attention. 
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